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When a beam of light is reflected by a smooth surface its behavior deviates from geometrical 
optics predictions. Such deviations are quantified by the so-called spatial and angular Goos-Hanchen 
(GH) and Imbert-Fedorov (IF) shifts of the reflected beam. These shifts depend upon the shape 
of the incident beam, its polarization and on the material composition of the reflecting surface. In 
this article we suggest a novel approach that allows one to unambiguously isolate the beam-shape 
dependent aspects of GH and IF shifts. We show that this separation is possible as a result of some 
universal features of shifted distribution functions which are presented and discussed. 



I. INTRODUCTION 

The reflection and refraction of light at the plane sur- 
face of a transparent medium, such as the still water of a 
placid lake, are well-known phenomena. As a matter of 
fact, the law of reflection was already enunciated by Eu- 
clid of Alexandria (fl. 300 BC) in his book Catoptrics^ 
Nowadays, extended expositions of this subject may be 
found in standard optics textbooksi^ 

Less well-known is the fact that when the light has the 
form of a beam with a finite transverse (i.e., orthogo- 
nal to the direction of propagation) extent, the reflected 
beam appears shifted with respect to a ray traveling par- 
allel to the beam but reflected according to the laws of 
geometrical optics, as opposed to the laws of wave optics 
obeyed by the finite beam. A cartoon-like representation 
of this phenomenon is presented in Fig. 1. 




FIG. 1. Illustration of a beam-shift process. The incident 
beam travels along with the ray represented by the unit vec- 
tor Zi. Reflection transforms Zi according to the laws of ge- 
ometrical optics into Zr = Zi — n (n ■ Zi), where n is a unit 
vector normal to the reflecting surface. The reflected beam 
travels parallel the unit vector fir which is both displaced and 
tilted with respect to Zr. Aqh and Aip denote, respectively, 
the projections onto the plane of the surface of the beam dis- 
placements parallel and orthogonal to the plane of incidence 
spanned by ii and fi. 



Such a shift may have both a spatial and an angular 
character, namely the refiected beam may be both dis- 
placed and deflected with respect to the geometrical ray. 



Depending on whether these spatial and angular shifts 
occur in either the plane of incidence or in a direction or- 
thogonal to this plane, they are denoted as either Goos- 
Hanchen (GH) or Imbert-Fedorov (IF) shifts. Analogous 
shifts also occur for the refracted beam. 

All these shifts are very small compared to the trans- 
verse extent of the beam, typically a few wavelengths the 
spatial shifts, and from micro- to milli-radians the angu- 
lar ones. Interestingly enough, the first who conjectured 
that the center of a beam of light should manifest a tiny 
spatial shift in the plane of incidence when refiected, was 
Isaac Newton in the 17th centuryi^ 

Detailed theoretical derivations of both GH and IF 
shifts may be found in the current literature for surfaces 
of different types (dielectric, metallic, multilayered, et 
cetera), and for light beams of vary shapes (Hermite- 
Gauss, Laguerre-Gauss, Bessel, et cetera). The inter- 
ested reader may consult Refs. 5-8 for early expositions 
and Refs. 9-15 for more accessible modern treatments of 
the subject. From these studies it emerges that, given a 
specific reflecting surface, beams of different shapes yield 
shifts of different extent. For example, for an air-to-glass 
interface, a fundamental Gaussian beam produces qual- 
itatively and quantitatively different GH and IF shifts 
with respect to a Laguerre-Gauss beam impinging upon 
the same surface. This could appear not so surprising 
because a Laguerre-Gauss beam carries an intrinsic or- 
bital angular momentum while a fundamental Gaussian 
beam does notJ^ Thus, one might conclude that orbital 
angular momentum (0AM) of light affects the physics of 
beam-shift phenomena. More generally, since such 0AM 
is determined solely by the spatial profile of the beam, 
one may wonder to what extent the beam shape deter- 
mines the values of both GH and IF shifts. 

Unfortunately, prevailing discussions in the current 
literature tend to obscure the distinction between the 
beam-shape dependent and independent aspects of 
beam-shift phenomena because conventional evaluations 
of GH and IF shifts rest on the direct calculation of 
the centroid of the intensity distribution of the reflected 
beam. Such a distribution depends, as a whole, upon 
both the shape of the beam and the material composi- 



tion of the surface. Thus, beam-shape dependent and in- 
dependent characteristics of shift phenomena deceptively 
appear to be intrinsicahy tied. 

In this article we present a novel way to describe the 
GH and IF shifts, with the aim to provide significantly 
new insights about these phenomena. To this end, a sim- 
ple treatment that renders a clear separation between 
the beam-shape dependent and independent features of 
beam-reflection processes is presented and developed. 
Such separation naturally leads to find some "universal" 
characteristics of the beam-shape dependent facets of the 
problem. 

The rest of this paper is organized as follows: In Sec. 
II the conventional theory of GH and IF shifts is reviewed 
and expressions for the beam-shape independent parts of 
these shifts are explicitly given. In Sec. Ill and IV a new 
treatment based on the study of the functional shape of 
the beam is introduced and illustrated. Finally, in Sec. 
V, consequences of the preceding results are discussed. 



II. GOOS-HANCHEN AND IMBERT-FEDOROV 
SHIFTS 

Consider a well collimated monochromatic beam of 
light of wavenumber k propagating in the direction z^. 
The geometric details are shown in Fig. 2 below. 
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FIG. 2. Geometric representation of the system. The plane 
of incidence coincides with the plane of the figure. The unit 
vector fi is normal to the interface between media 1 (top) and 
2 (bottom). The unit vectors yi = y^ = y (not shown in 
the figure), are perpendicular to the plane of incidence and 
directed towards the reader. The Snell law of reflection sets 
6i = 9r = 9. The index of refraction n\ of the flrst medium is 
conventionally chosen equal to 1, and the index of refraction 
n2 of the second medium is n. 



The position vector r has Cartesian co- 
ordinates r = [x' Xi + y' Yi + z' ii)/k and 
T = {x'kr -\- yfr ~\- z2,r)/k with rcspcct to the refer- 
ence frames {xi,yi,Zi} and {kr,Yr,ir} attached to 
the incident and to the reflected beam, respectively. 
Hereafter all the distances x',y',... are dimensionless. 



The corresponding dimensional quantities are simply 
equal to x'/fc, y' /k^ et cetera. The real electric field 
associated with such a beam may be written as 

E'(a;',2/',/,i) = Re [A'(x', y', z') exp (-i/cci)] , (1) 

where c denotes the speed of light in vacuum and 

A'{x', y', z') = [a || x^ + ai_ y,) f{x', y' , z'), (2) 

is the so-called complex analytic signalr^ associated to 
the real signal E'(a;', y', z', t). In Eq. ^ and hereafter, 
subscripts "||" and "_L" refers, respectively, to the di- 
rections parallel and perpendicular to the plane of inci- 
dence. In optics textbooks such directions are often de- 
noted with the letters "p" and "s", respectively. The 
two complex amplitudes a ii and a ± define the polar- 
ization of the beam and the complex scalar function 
f{x', y' , z') fixes the spatial shape of the beam. For exam- 
ple, f{x',y',z') may either be a Hermite-Gauss function 
or a Laguerre-Gauss function, both types being widely 
utilized in optics.~ 

When the beam is reflected its analytic signal 
A'(a;', J/', z') changes into A(x, y, z), where^^ 

A{x,y,z) = a[|r||(6l)/(-.T + X||,y-r||,z)xr 

+ a^r^[e)f{-x + X^,y^Y^,z)yr, (3) 

and ?'||(6'), r±{0) are the Fresnel reflection coefficients 
evaluated at the incident angle Oi = 6. The four di- 
mensionless complex shifts X||,y|| and X^^Y±_ are given 
by.^ 
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These shifts depend upon the polarization of the incident 
beam via the complex amplitudes a|| and a±, and upon 
the characteristics of the surface via the (either real or 
complex) Fresnel coefficients r|| and r±. However, they 
are com,pletely independent from the shape of the beam 
f{x',y',z'). The reader should notice and keep in mind 
this important fact. Equations ^ and @ are valid un- 
der the assumption that the reflection does not strongly 
distort the beam. Such deformation occurs, for exam- 
ple, when the beam impinges at the Brewster angled for 
which one has rn — and Yu becomes ill-defined; in this 
case a special treatment is requiredi^^ 

For monochromatic light the intensity of the reflected 
beam is given byji 
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where Eq and /xq are the vacuum permittivity and 
permeabihty,— respectively, and 

I 1 2 I 1 2 

-^11 = rinil \fi-x + ^hy-~Yn,z)\ , (6a) 

I 1 2 I 1 2 

I± = \a±rj_\ \f{-x + Xj_,y-Y^,z)\ , (6b) 

are the intensity distribution functions for the reflected 
light with polarization parallel and perpendicular to the 
plane of incidence, respectively. 

The center of the reflected beam may be evaluated as 
the centroid in the xy plane of the total intensity distri- 
bution /|| +/^ . The position of such centroid with respect 
to the origin x = 0, y = on the plane z = const., de- 
fines both the spatial and angular GH and IF shifts in the 
following way. Given the Cartesian components xi = x 
and X2 = y oi the centroid of the beam defined as: 



{xa){z) 



Xa {l\\ 



dx dy 



(a = 1,2), (7) 



[h+Ii 



dxdy 



then the spatial (A) and angular (0) GH and IF shifts 
are conventionally defined as: 



Agh = (x)(0), 
AiF = (2/)(0), 



d{x){z) 
fcH — — , (oa) 
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(8b) 



Note that in Eq. (O and hereafter the integration sym- 
bol "/ " always denotes definite double integrals over the 
whole xy plane with z kept constant. Moreover, the fact 
that 0GH and 0if are strictly independent from z (see, 
e. g., Ref. 22) is also worth noticing. 

It is easy to see that after a straightforward algebraic 
manipulation Eq. ([7]) can be recast in the following sim- 
ple form: 

{Xa){z) = 'W\\{Xa)\\{z) + W^{Xa) x_{z), (9) 

where the weight coefficient wn is defined as: 
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and the relative shift as: 



{xc)\\{z) 



Xa\f{-x + X\\,y -Y\\,z)\ dxdy 

I 1 2 

\}{-x + X\\,y-Y\\,z)\ dxdy 



(11) 



The corresponding formulae for w^^ and {xa) ^{z) can 
be obtained from Eqs. (|10j|lip by exchanging the sub- 
scripts "II" and "_L" everywhere. Equation (|10p has been 



obtained by exploiting the fact that since the double in- 
tegration extends over the entire xy plane, one has: 

\f{-x + X\\,y-Y\\,z)\ dxdy = 

J\fi'X + X^,y-Y^,z)\^dxdy. (12) 

Equations (jlOp and dill) have a simple physical interpre- 
tation. The non-negative weight coefficients wh and w± 
give the fraction of the reflected light which is linearly 
polarized in the directions parallel and perpendicular to 
the plane of incidence, respectively. The relative shifts 
{xa)\\{z) and {xa)±{z) are the first moments of the dis- 
tributions /|| and /_L, and represent the displacements 
of the intensity distributions that would be detected be- 
hind a polarizer oriented along x^ and y^, respectively, 
at distance z from the origin. 

At this point, conventional calculations require substi- 
tution of Eqs. (HI into Eq. pT|) and subsequent straight- 
forward integration, either analytical for a few simple 
cases, or numerical for the more complicated ones.— 
However, as announced in the Introduction, such pro- 
cedure leads to a confusion, in determining GH and IF 
shifts, between the role played by the functional shape 
f{x, y, z) of the beam on one side, and the role played by 
the beam-shape independent shifts ^||,i^|| and X^,Y±^ 
on the other side. This issue will be clarified in Sec. HI 
and practically illustrated in Sec. IV. 

III. FUNCTIONAL SHIFTS 

In order to elucidate the problem posed at the end 
of the previous section, a seemingly more abstract issue 
is illustrated here to the reader in some detail: Given 
an arbitrary complex function f{x, y, z) square integrable 
over the whole xy plane, calculate the first moments {x) 
and (y) of the shifted distribution 



F{x,y,z) = \J{x^X,y-Y, 



(13) 



where X,Y £ <C and \X\, \Y\ <^ 1.^* In other words one 
is required to calculate the following quantities: 



{x) = 



(y) 



xF(x, y, z) dxdy 
F{x,y, z) dxdy 

yF(x, y, z) dxdy 
F{x,y, z) dxdy 



(14a) 



(14b) 



The complex shifts X and Y of the argument of the func- 
tion are written as: 



X = Xi+iX2, 
Y^Y^+iY2, 



(15a) 
(15b) 



where X„,F„ G R, with n G {1,2}. A first-order Taylor 
expansion about X = = Y furnishes: 

F{x,y,z) = |/(x,2/,z)p 
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where hereafter / is a shorthand for f{x, y, z), not to be 
confused with f{x—X, y—Y, z). As usual, "Re" and "Im" 
denotes, respectively, the real and the imaginary part of 
a complex number: Re(a + ib) = a and Im(a + ib) — b, 
with a, 6 S R. Without loss of generality, we assume that 
the function / drops to zero at infinity and that both 
the X- and y- moment of the distribution j/p are zero, 
namely: 



|/(±c^, y,z)f = 0= \fix, ±c^, z)|2, (17a) 

x\f{x,y,z)fdxdy = 0= y\f{x,y,z)fdxdy. (17b) 



By using integration by parts, it is not difficult to show 
that, up to first order terms in Xn, Yn, one obtains after 
substitution of Eq. ([H]) into Eqs. ([Ti)) : 



(x) = Xi + 2 Im 



(y) ^ Fi + 2 Im 
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(18a) 



(18b) 



One should appreciate the fact that these formulae are 
perfectly general, as they hold their validity for any 
square-integrable function f{x,y,z) fulfilling conditions 

The meaning of Eqs. (fT8|) becomes crystal clear when 
they are put in the following suggestive matrix form: 



(19) 



where the elements aap of the matrix A = [aq^], with 
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an ai2 

^21 022 




' X2' 
.^2 _ 



a,/3 — 1,2, have been defined as: 

Im f* [xa T. — I / dxdy 

aai3 = 2 j: . (20) 

Iff dxdy 



The first term in the right side of Eq. ([T^ simply states a 
trivial fact: If one shifts a distribution, initially centered 
in the origin of the xy plane, by the real amounts Xi 
and Yi in the x and y directions, respectively, then the 
shifted distribution will be centered in P = (Xi, Yi). 

Much more interesting is the second term. The diago- 
nal elements of the 2x2 matrix A connect the imaginary 
parts of the complex shifts X and Y with the real dis- 
placements of the distribution. Since from Eq. ([^0]) it 
trivially follows that A must be a function of z solely, 
this means that X2 and I2 rnust somehow describe the 
propagation properties of the shifts (x) and (y) , respec- 
tively. It will be explicitly shown in Sec. IV that this is 
indeed the case for functions / describing paraxial beams 
of light, and that an and 022 are responsible for the an- 
gular GH and IF shifts, respectively. 

The non-diagonal matrix element ai2 (021) mixes the 
real part of the shift X (Y) in the x (y) direction with 
the imaginary part of the shift Y (X) in the y (x) direc- 
tion. The existence of these nonzero off-diagonal matrix 
elements is a remarkable universal feature of shifted dis- 
tributions. In fact, although the exact amount of mixing, 
namely the precise values of 012 and 021, depends on the 
functional shape of /, its existence is largely independent 
from the form of /. This seemingly purely mathemati- 
cal property of distributions has actually profound con- 
sequences on the physics of beam-shift phenomena. In 
fact, such a mixing has been actually predicte d^^'^^ and 
experimentally observed^S for vortex beams. Thus, for a 
better understanding of this issue, in the following Sec. 
IV the perfectly general mathematical analysis of func- 
tional shifts given above will be specialized to the physi- 
cal case of vortex beams, namely beam carrying 0AM. 



IV. EXAMPLE: FUNCTIONAL SHIFTS OF 
VORTEX BEAMS 

The goal of this section is to illustrate the use of the 
formula ([T9| by calculating the coordinates of the cen- 
troid of a reflected vortex beam. According to Eq. (jH]) 
the shifted beam is described by one or two functions of 
the form /(— x + X,y — Y, z), with 



fix,y,z) 



^imcj) 



9{r,z), 



(21) 



where to is a real number, not necessarily an integer one, 
and X = r cos 4>, y — r sin cf) are the polar coordinates in 
the xy plane. Equation (j21[) . with integer to, furnishes 
a suitable representation of optical beams with m units 
of 0AM along the direction of propagation z. For ex- 
ample, the well-known Laguerre-Gauss paraxial modes 
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tional form 



u'^((f>, r, z) of the electromagnetic field, have the func- 



9ir,z) 



B with 
2pl 



TT {p+ \m\)l 



1/2 



exp 



X exp 
1 



-z (2p 



w{z) p 



2z-iL 
1) arctan 



2r' 



w'^{z) J \ w{z) 



(l 

|rn| 



(22) 



where wq is the beam waist, L = k'^wff/2 is the Rayleigh 



range, w{z) = fcwoi/l + z"^ j LP- is the running beam 

waist, and L^ {x) denotes the generalized Laguerre 
polynomial^i of order {ni,p), with p = 0,1,2,... and 

m = 0, ±1, ±2, The first few Laguerre-Gauss modes 

intensity patterns are pictured in Fig. 3 below. 
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FIG. 3. Some low-order Laguerre-Gauss modes intensity pat- 
terns \ul;^{(ji,r,z)\ . 

For the sake of simplicity, hereafter X, Y denote either 
X|| , Yji or X± ,Y±, we leave it unspecified. Correspond- 
ingly, for example, (x) denotes either (a;)|| or (a;)_L. Thus, 
from Eqs. (|20 l) - (l2T|) . we can obtain the matrix A by cal- 
culating the following integrals: 



\fi-x,y,z)\'^dxdy = 27r / \g{r, z)\'^rdr 



(23) 



and 



/ ■^* \'dx) ■^^^^^ "" ^ / 9*{r,z) {r—j g{r,z)rdr 



N, 



/* f a;^ j / dxdy = - inm / \g{r,z)\'^ rdr 
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(24a) 



(24b) 



where / now is a shorthand for f{—x,y,z). The last 
remaining two integrals are not independent from the 
previous ones, being: 



/* 



2/g^ ) fdxdy 



N. 



r[y—]fdxdy= -iM. 



Finally, substitution of Eqs. dM]) into Eq. ([T 
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(25a) 
(25b) 
yields: 

(26) 



where 



C(^) 



Im / g*{r,z)[r—]g{r,z)rdr 



\g{r,z)\'^rdr 



(27) 



determines how the shifted distribution "propagates" 
along z. In physical terms, the real function C,{z) gen- 
erates the angular GH and IF shifts. This may be easily 
seen by deriving both sides of Eq. (PB)) with respect to 
z and comparing such result with the right sides of Eqs. 



(28) 



Of course, the precise form of C,{z) can be given only 
when g(r, z) is precisely specified. Thus, for example, for 
a Laguerre-Gauss paraxial beam g(r, z) of the form 
one obtains from Eq. (P7| : 
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where, for the sake of clarity, the value p = has been 
chosen. When this expression is substituted into Eq. (^5)) 
the following result is obtained: 



(x) 9iXi-mY2 + - (1 



A)X2, 



{y)=Yi+mX2 + -{l + \m\)Y2. 



(30a) 
(30b) 



If TO = the Laguerre-Gauss beam reduces to a fun- 
damental Gaussian beam, and Eqs. (I30p become: 



(x)" - Xi + -X2, 

(2/>° ^yi + 2Y2, 



(31a) 
(31b) 



where hereafter the superscript "0" denotes quantities 
relative to a fundamental Gaussian beam. This leads to 
the unambiguous identification 



Xi 


- A° 


X2 = ieg;H, 


(32a) 


Yi 


- AO 

— ^IF' 


Y2 = LQ%, 


(32b) 



where Eqs. ([8]) have been used. Finahy, by using Eqs. 
([5^ into Eqs. ([50]) . one easily obtains: 






(y) = A?F 



n(Le?F)+z(l + |77i|)eg;H, (33a) 
.{Le%^)+z{l + \m\)e%. (33b) 



This outcome coincides, except for a few minor discrep- 
ancies due to the use of a different notation, with the 
well-known results for vortex beamsi ^°i^^'^^ From Eqs. 
((33)l it appears that 0AM of light has a twofold effect 
upon GH and IF shifts. The first one is an algebraic 
mixing proportional to m between spatial GH (IF) and 
angular IF (GH). The second effect is an increase of the 
angular shift by the amount \m\. As revealed by the ap- 
proach presented in this article, these features are univer- 
sal for Laguerre-Gauss beams, namely they are indepen- 
dent upon both the beam polarization and the material 
composition of the reflecting surface. 

Thus, Eqs. (|33l) and ((T^ answer the two questions 
posed in the introduction, namely: a) how the 0AM of 
light affects GH and IF shifts, and: b) to what extent 
these shifts are determined by the shape of the beam 
(see discussion after Eq. (pUp. 



V. CONCLUDING REMARKS 

The main result presented in this paper perhaps con- 
sists in the following recipe for the calculation of GH and 



IF shifts: Use the conventional theory to calculate the 
beam-shape independent complex beam shifts X|| , 1]| and 
X± , Y± , for your peculiar surface. They depend upon 
both the polarization of the beam and the material com- 
position of the reflecting surface, but are totally inde- 
pendent from the shape of the beam. Then, by using 
Eq. ((20|). calculate the matrix A for your speciflc beam. 
Finally, combine the previous results into Eq. (J19p to 
obtain the desired shifts. 

Apart from its great practical usefulness, this new ap- 
proach has the important conceptual advantage over the 
traditional beam-shift computational methods, that it 
permits to isolate in a completely unambiguous manner 
the beam-shape dependent contributions to the shifts. 
This separation is relevant whenever one wants to under- 
stand how the degrees of freedom carried by the wave- 
front of the beam (as, e.g., the 0AM of the beam), enter 
in the numerous beam-shift physical phenomena. 
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